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Abstract 

We review various attempts to localize the discrete spectra of semirelativistic Hamiltonians 
of the form H = f3\/m 2 + p 2 + V(r) (defined, without loss of generality but for definiteness, 
in three spatial dimensions) as entering, for instance, in the spinless Salpeter equation; every 
Hamiltonian in this class of operators consists of the relativistic kinetic energy f3\Jm 2 + p 2 , 
where (3 > allows for the possibility of more than one particles of mass m, and a spherically 
symmetric attractive potential V(r), r = |x|. In general, accurate eigenvalues of a nonlocal 
Hamiltonian operator can only be found by the use of a numerical approximation procedure. 
Our main emphasis, however, is on the derivation of rigorous semi-analytical expressions for 
both upper and lower bounds to the energy levels of such operators. We compare the bounds 
obtained within different approaches and present relationships existing between the bounds. 
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1 Introduction: The "spinless-Salpeter" Hamiltonian 

We investigate the discrete spectrum — that is, the set of eigenvalues E n £ below the onset of 
the essential spectrum, corresponding to bound states characterized by the radial quantum 
number n — 1, 2, 3, . . . and the angular-momentum quantum number t = 0, 1, 2, . . . — of the 
semirelativistic "spinless-Salpeter" Hamiltonian for particles of mass m and momentum p, 

H = ^m 2 + p 2 + V{r) , r=|x|, (3 > , (1) 

where V(r) is an attractive central potential in three spatial dimensions. As is evident from 
the presence of the relativistic kinetic-energy operator fi^pmF^Pp 1 (with (3 > allowing, 
e.g., for more than one particle), this Hamiltonian may be regarded as the generalization of 
the nonrelativistic Schrodinger Hamiltonian towards relativistic kinematics. The eigenvalue 
equation of the operator (1), called the "spinless Salpeter equation," arises as a well-defined 
approximation to the Bethe-Salpeter formalism for the description of bound states within a 
relativistic quantum field theory. Our main goal is to discuss bounds on the energy levels E. 

2 Envelope theory for energy bounds 

The envelope theory, developed in Refs. [1,2], constructs bounds on the discrete eigenvalues 
E n £ of a given operator H by comparing the spectrum of H with the spectrum of a suitably 
formulated "tangential problem," for which some spectral information is known (or may be 
obtained with considerable more ease than for the original problem). More specifically, when 
applied to the (semirelativistic) Hamiltonian H defined in Eq. (1) [3, 4], the envelope theory 
compares H with a corresponding "tangential Hamiltonian" 

H = j3^Jm 2 + p 2 + vh(r) , v > , 

involving some "basis potential" h(r). Under the assumption that the interaction potential 
under consideration, V(r), is a smooth transformation V(r) = g(h(r)) of the basis potential 
h(r), where the function g(h) has definite convexity, this comparison yields rigorous bounds 
on the discrete eigenvalues E n £ of H. More precisely, if the function g(h) is convex, that is, 
g"{h) > 0, we obtain lower energy bounds; if the function g(h) is concave, that is, g"(h) < 0, 
we obtain upper energy bounds. If the eigenvalues E of H are not known (analytically), the 
envelope theory may take advantage of the knowledge of suitable bounds on the eigenvalues, 
namely, of lower bounds on E n £ for convex g{h) or of upper bounds on E n £ for concave g(h). 

Suppressing the quantum numbers n, £ identifying the bound state under consideration, 
we denote, for a given value of the coupling v entering into the "tangential Hamiltonian" H, 
the eigenvalues of H or the appropriate — in the above sense — bounds to the latter by e(v). 
The strict bounds on the discrete eigenvalues E of H predicted within the framework of the 
envelope technique may then be summarized in form of the "principal envelope formula" [4] 

E » mm[e(v) - ve'(v) + g{e\v))\ , (2) 

where the sign of approximate equality is used to recall that, for a definite convexity of g(h), 
the envelope theory yields lower bounds for convex g{h) and upper bounds for concave g(h). 

For the tangential problem posed by the (semirelativistic) "tangential Hamiltonian" H, 
spectral information is, at present, available if the basis potential h(r) is either the Coulomb 
potential, that is, h(r) = — 1/r, or the harmonic-oscillator potential, that is, h(r) = r 2 . The 
energy bounds resulting in these two cases are discussed, in turn, in Sec. 4 and Sec. 5. The 
case of the interaction potential V(r) in Eq. (1) being the harmonic-oscillator potential may 
be traced back to a nonrelativistic Schrodinger problem; this is studied for its own in Sec. 3. 
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3 Harmonic-oscillator potential 

The "relativistic harmonic-oscillator problem," defined by the special case of V(r) being the 
harmonic-oscillator potential V(r) = vr 2 , may be analyzed by observing [5, 6, 7, 8, 9] that in 
momentum-space representation the semirelativistic Hamiltonian H = fi^Jrri 1 + p' 2 + vr 2 of 
Eq. (1) simplifies to the (nonrelativistic) Schrodinger operator H given by 

H = -vA + V(r) , V(r) = [3Vm 2 + r 2 . (3) 

Thus the eigenvalue equation for H reduces to an easier-to-treat nonrelativistic Schrodinger 
problem with an effective interaction potential V(r) which is reminiscent of the square root. 



3.1 Upper bounds on harmonic-oscillator energy levels 

The potential V is a concave transform of a harmonic-oscillator potential. Thus we may find 
upper bounds on the eigenvalues E ni (v) of H (cf. Eq. (2.2) of Ref. [3] or Eq. (12) of Ref. [4]): 



E n i(v) < min 

r>0 



P nl {2) = 2n + i- 1 -. (4) 



The (dimensionless) numbers P n e(2) introduced here are related to the eigenvalues S n £ of the 
Hamiltonian H = p 2 + vr 2 by E n t = 2y/vP n c(2). 

Upper bounds on the eigenvalues of self-adjoint operators bounded from below may also 
be found by combining the minimum-maximum principle [10, 11] with appropriate operator 
inequalities. (For details, see, e.g., Sec. 3.4 of Ref. [5], Sec. 3.5 of Ref. [6], Sec. 2.4 of Ref. [7], 
or Appendix A of Ref. [12].) Exploiting the obvious positivity of {y/rri 1 + p 2 — /i) 2 , where \x is 
an arbitrary real parameter (with the dimension of mass), we obtain a set of inequalities for 
the square- root operator ^m 2 + p 2 entering into the kinetic energy of the Hamiltonian (1); 
cf. also Ref. [13]: 

yf^T^< p2 + f + /i2 V„>0. 

_ LL 

(This operator inequality may also be found by constructing the tangent line a(Pi)p 2 +fr(p 2 ) 
to the square root \frri 1 + p 2 at the point of contact p 2 = /i 2 — m 2 ; see Sec. 4 of Ref. [4].) 
Amending this result by the interaction potential V(r), one finds that the spinless-Salpeter 
Hamiltonian (1) satisfies 

H < ff NR (/i) = (3 p2 + ^ + ^ + V(r) V/i>0. 

The min-max principle then tells us that the discrete eigenvalues E n £ of the semirelativistic 
Hamiltonian H of Eq. (1) are bounded from above by the corresponding discrete eigenvalues 
E-^R^ntifJ') of the nonrelativistic Hamiltonian H^n(n) and consequently also by the minimum 
E n £ with respect to /i of all these upper bounds [5]: 

Efd < E nl ee mm E NR n £(/_i) . 

Applying these general findings to the case of the harmonic-oscillator potential, one ends up 
with the following upper bounds for all energy levels of the relativistic harmonic-oscillator 
problem (see Appendix A of Ref. [9]; for the ease of comparison, we express here the bounds 
of Ref. [9] in terms of P n t{2) rather than in terms of E NRjn e(/j,)): 



E n e(v) < E nl ee mm 




? (2) = 2n + £-i. (5) 
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Despite their seemingly different form, the upper bounds in Eqs. (4) and (5) are, in fact, 
identical, in the sense that the two functions on the right-hand side of these inequalities have 
their minimum at the same critical point and take the same miminal value at this point. 

In order to demonstrate this equivalence, let us start from the "min-max" expression (5) 
for these upper bounds. For simplicity of notation, we suppress for the moment the quantum 
numbers n, £ which identify the bound states under consideration. The upper bound E then 
reads 

m 2 + fi 2 [J3v c 



E = min 

M>0 



— + J—2P 

1 1\i V 2// 



We introduce the abbreviation q = p 2 , define the function f[[q) = \/m 2 + q, which satisfies 

/'(?) = 2^y, f 2 (q)=m 2 + q, 

and identify the (positive but otherwise arbitrary) parameter \i according to /i — f(q). Next, 
we introduce a function F by the definition F(v) = 1P\Jv. A couple of algebraic steps then 
recasts the upper bound E into the form 

E = mm{F(0vf') + 0(f-qf')] . 

The critical points of this latter expression are obtained by differentiation with respect to q. 
Since f"{q) ^ 0, this yields 

q = vF'(0vf'(q)) 

and 

E = F(0vf(q)) - 0vf{q)F'{0vf{q)) + 0f{q) . 

Let us define a new problem, related to the above one by a Legendre transformation [14], 
which has the same minimum and the same critical point, but the function to be minimized 
is different. Let u(q) = 0vf'(q) > and consider the energy expression 

E = min[F(ti) - uF\u) + 0f(vF'(u))] . 

Since F"{u) ^ 0, we find that the critical point is given by u = 0vf'(vF'(u)). But, given the 
definition of u, this equation is equivalent to the earlier critical equation q = vF'(0vf'(q)). 
Meanwhile, the critical value E is also the same as the critical value E because we have 

E = F(0vf(q)) - 0vf{q)F'{0vf{q)) + 0f(q) = E . 

Having shown that E yields the same results as E, we now transform E into the form (4). 
First of all, we note that F(u) —uF'{u) = Py/u. Next, we define a new dummy minimization 
variable r > by 



2 

vr 



From this, it follows that 

P P 2 
F\u) = — = — . 
\Ju vr z 



Finally, remembering the definition f(q) = \Jm 2 + q of f(q), we may transform E into the 
expression on the right-hand side of the inequality (4); this establishes the equivalence of the 
upper bounds (4) and (5). 
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3.2 Lower bounds on harmonic-oscillator energy levels 



The potential V in Eq. (3) is a convex transform of a linear potential. Therefore we may find 
lower bounds on the eigenvalues E ni (v) of H (cf. Eq. (2.2) of Ref. [3] or Eq. (12) of Ref. [4]): 



E n Av) > min 



(3\ m 2 + 



(6) 



where the (dimensionless) numbers P n (,{l) introduced here are related to the eigenvalues £ n i 
of the Hamiltonian H = p 2 + v r by 



£ni — 3 



2/3 



for the lowest states {n = 1, 2, . . . , 5, £ = 0, 1, . . . , 5) these numbers may be found in Table 1. 



Table 1: Numerical values of the energy-related numbers P n i(l) for the lowest energy levels. 



n 


i 




1 





1.37608 


2 





3.18131 


3 





4.99255 


4 





6.80514 


5 





8.61823 


1 


1 


2.37192 


2 


1 


4.15501 


3 


1 


5.95300 


4 


1 


7.75701 


5 


1 


9.56408 



n 


i 


Pni(l) 


1 


2 


3.37018 


2 


2 


5.14135 


3 


2 


6.92911 


4 


2 


8.72515 


5 


2 


10.52596 


1 


3 


4.36923 


2 


3 


6.13298 


3 


3 


7.91304 


4 


3 


9.70236 


5 


3 


11.49748 



n 




P n e(l) 


1 


4 


5.36863 


2 


4 


7.12732 


3 


4 


8.90148 


4 


4 


10.68521 


5 


4 


12.47532 


1 


5 


6.36822 


2 


5 


8.12324 


3 


5 


9.89276 


4 


5 


11.67183 


5 


5 


13.45756 



Lower bounds on the eigenvalues of (self-adjoint) Hamiltonians involving the relativistic 
kinetic energy y/m 2 + p 2 may also be found by taking into account the obvious positivity of 
[p 2 £ 2 ~ rn 2 (l — £ 2 )] 2 for an arbitrary real parameter £, which may be converted into the set 
of operator inequalities 

yW^ + p 2 ^ Ipl^/l-^ + m^, 0<C<1. (7) 

Applying these general findings to the semirelativistic Hamiltonian H of Eq. (1), with V(r) 
being the harmonic-oscillator potential, we obtain lower bounds on all discrete energy levels 
E n £ of the relativistic harmonic-oscillator problem (for details, see Appendix A of Ref. [9]): 

E ne (v) > E nl = max[/3£m + (3 2 / 3 (l - £ 2 ) 1/3 ^; m = 0)] , (8) 

where E n g{y\ m — 0) denote the corresponding discrete eigenvalues of the operator |p| +vr 2 , 
i.e., of the Hamiltonian (1) for (3 — 1 and a vanishing particle mass m. For states with t = 0, 
for instance, these "zero-mass" energy levels are given, in terms of the zeros z n of the Airy 
function M(z) [15] {-z n = 2.33810, 4.08794, 5.52055, . . .), by E n0 (v;m = 0) = -v l l z z n [9]. 
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As in Sec. 3.1 we observe that the lower bounds in Eqs. (6) and (8) appear, at first sight, 
to be rather different. A closer inspection, however, reveals that these bounds are identical. 

In order to present an analytical proof of the equivalence of the lower bounds (6) and (8) 
on the energy eigenvalues of the relativistic harmonic-oscillator problem, we recall that the 
envelope lower bounds in (6) are based (in the momentum-space representation of H) on the 
observation that, since its second derivative with respect to |p| is positive, the square-root 
operator y/m 2 + p 2 of the relativistic kinetic energy is bounded from below by tangent lines 
a(t)+b(t)\p\] this may be expressed by the geometrical inequality ypmF^Pp 1 > a.(t)-h£»(t)|p>|, 
where, by elementary calculus, 

TV? t 

a(t) = - , b(t) = - , 

V in 2 + t 2 v m 2 + t 2 

and |p| = t is the point of contact between the square root of the relativistic kinetic energy 
and its straight-line lower approximation. The equivalence of the above two expressions for 
the harmonic-oscillator lower bound is then shown by a simple change of variable from t to £ 
given explicitly by 

t _ <t) 



which, after a little algebra, leads to b(t) = \/l — £ 2 . For these coefficient functions a(t) and 
b{t) the above geometrical inequality becomes exactly the lower bound (7) on the relativistic 
kinetic energy. Consequently, with the basic inequalities being identical, the resulting lower 
energy bounds have to be identical too. 



4 Convex transform of Coulomb potential 



For the "relativistic Coulomb problem" posed by the spinless-Salpeter Hamiltonian (1) with 
V(r) being the Coulomb potential V(r) = —v/r, for coupling constants v smaller than (3v c , 
that is, v < (3v c , where the critical value v c of the Coulomb coupling constant v is given by 

2 



7T 



a lower bound to the ground-state energy eigenvalue E of H (the bottom of the spectrum 
of H) has been derived by Herbst [16]: 



£ > Pm 



1 - 



av 
J, 



a = 



(9) 



For some part of the allowed range of the coupling v, this lower bound has been improved by 
Martin and Roy [17]: 



E > (3m 



\ 



i + 



\ 



V< 2- 



(10) 



In Ref. [4], we recast this Martin-Roy bound into a class of lower bounds on the spectrum of 
H all of which are of the form of Herbst's bound (9) [but, of course, slightly weaker than the 
Martin-Roy bound (10)]; the new lower bounds are parametrized by a real parameter a > 1: 



£ > Pm 



6 



The constraint on the Coulomb coupling v in Eq. (11) automatically guarantees v < j3/cr; 
this latter constraint is a necessary and sufficient condition for the reality of the lower energy 
bounds in Eq. (11). 

From the lower bounds (9) or (11) on the spectrum of the relativistic Coulomb problem, 
for the class of potentials V(r) which are convex transforms V(r) = g(h(r)), g" > 0, of the 
Coulomb potential h(r) = — l/r an ("absolute") envelope lower bound on the ground-state 
eigenvalue Eq (and, hence, on the entire spectrum) of the spinless-Salpeter Hamiltonian (1) 
may be found (cf. Eq. (16) of Ref. [4]): 



Eq > min 

r>0 



pJm 2 + ^ + V(r) 



, v < (3vp 



(12) 



Here we introduced a parameter P by defining P = 1/er. Accordingly, the boundary value Vp 
of the Coulomb coupling constant v is given, when arising from demanding the Hamiltonian 
(1) to be bounded from below, by the critical coupling v c , i.e., vp = v c , or, when arising from 
the region of validity of our new family of Coulomb lower bounds (11), by 



v P = PVl - P 2 < - ■ 



(13) 



It is easy to convince oneself that the Coulomb lower energy bound on the right-hand side of 
the inequality (12) is a monotone increasing function E_(P) of the parameter P. To this end, 
we rewrite the inequality (12) in the form 



Eq > E(P) = min F(r,P) , 

r>0 



with 



F{r,P) = (3Jm 2 + ^ + V{r) . 



Remembering that V(r) is the transform V(r) = g(h(r)) of h(r) = —l/r, and assuming that 
g'(h) > 0, the critical point f(P) (corresponding to the minimum of F(r, P) for a given P) 
may be found, from 

dF(r, P) 



by solving 



dr 
[3P 2 



= 0, 



r=f(P) 



= g'(Hf)) . 



\Jm 2 f 2 + P 2 

The substitution of the solution f(P) of the latter equation in F(r, P) gives the lower bound 
E_(P) = F(f(P), P); here we are interested in the total derivative of E_(P) with respect to P. 
Formally, this derivative would be given by 



dE(P) _ dF(f(P),P) dF{r,P) 



dP 



dP 



dr 



df(P) 
dP 



r=f(P) 

At the critical point, however, the first factor of the second term vanishes and one is left with 

dE{P) dF(f(P),P) pp 



dP 



dP 



m 2 f 2 



+ P 2 



> . 



From this we conclude that the lower bound E(P) is a monotone increasing function of P. 
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The existence of the upper bound v < (3v P < (3/2 on the Coulomb coupling constant v is 
reflected in a corresponding constraint on the parameters of the interaction potential V(r). 
The envelope theory constructs a "tangential potential" to V(r) at a point of contact r = t. 
For a potential V(r) which is a convex transform V(r) = g(h(r)) of the basis potential h(r), 
this tangential potential forms a lower bound to V(r); for /i(r) being the Coulomb potential, 
h(r) = —l/r, this lower bound is a Coulomb potential with the effective coupling constant 
v {t) = t 2 V'(t), shifted by the constant V(t) + tV'(t): 

t 2 V'(t) 

V(r) > V(r) = + V(t) + tV'(t) . 

The relativistic Coulomb problem is posed by the "Coulombic" Hamiltonian 

H c (v) = (3\/m 2 + p 2 - - , v < (3v c , 

V r 

which is nothing else but the tangential Hamiltonian H with h(r) = —l/r. In the notation of 
Sec. 2, the eigenvalues Eq(v) of Hq(v) are bounded from below according to Eq(v) > e(v), 
with e(v) given by Herbst's bound (9), our new bound (11), or the Martin-Roy bound (10). 
Expressed in terms of the Coulombic Hamiltonian H c (v), the semirelativistic Hamiltonian 
(1) thus satisfies the inequality H > Hc(v(t)) + V(t)+tV (t); its eigenvalues E are bounded 
from below according to 

E > E c , (v(t)) + V(t) + tV'(t) > e(v(t)) + V(t) + tV'(t) . 

The optimized lower bound on the discrete spectrum of H is then found by maximizing the 
expression on the right-hand side of this latter inequality with respect to the point of contact 
r — t. The equation which determines the critical point t has to be solved by observing that, 
at the critical point r = i, the effective Coulomb coupling v(t) still has to satisfy v(t) < (3vp; 
this nontrivial requirement imposes the announced "Coulomb coupling constant constraint" 
on the coupling parameters entering in the potential V(r). 

It goes without saying that, when applying our "Herbst-based" lower energy bounds of 
Eq. (11), care has to be taken in order to assure that the chosen values of the parameters in 
the potential V(r) under consideration respect this Coulomb coupling constant constraint. 

Let us illustrate these general considerations by applying them to a Coulomb-plus-linear 
or (in view of its shape) "funnel" potential 

V(r) = -- + c 2 r , ci>0, c 2 > . (14) 
r 

In this case, the effective Coulomb coupling v(t) explicitly reads v(t) = ci + c 2 t 2 and, for the 
set of "Herbst-like" lower energy bounds e(v) given by Eqs. (9) and (11), the critical point t 
is, after some algebra, determined by the relation 

v(t) = = . 
\]m 2 t 2 + P 2 

The above Coulomb coupling constant constraint becomes an upper bound on a particular 
linear combination of the two involved coupling parameters c\ and c 2 (cf. Sec. 6 of Ref. [4]): 



p2 ,p2 \ 
C l H 9 -5" _ 1 C2 < f jVp 

m z \vf, J 



(15) 
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The Coulomb lower bounds (12) on the spectrum of the semirelativistic Hamiltonian (1) 
may be optimized by seeking that value of the parameter P = 1/a which, for given values of 
the mass m of the bound-state constituent (s) and of the coupling parameters c\ and C2 in the 
Coulomb-plus-linear potential (14), solves the Coulomb coupling constant constraint (15). 
Taking into account the definition (13) of the boundary value vp of the Coulomb coupling v, 
this optimum value of P is obtained as the solution of the relation 



C2 sin t 



cos 2 1(/3 sin t cos t — c\) 



m 



P = sint 



(16) 



The set (12) of envelope lower bounds distinguished from each other by the parameter P 
may be improved somewhat by taking into account the Martin-Roy lower bound (10) on the 
energy spectrum of the relativistic Coulomb problem. To this end, let us denote the energy 
bound on the right-hand side of the inequality (10) by e(v): 



e(v) = f3m 
Performing the change of variables 



1 






\ 2 


1 +N 


-(Si 



P IT 

v = £ sin 2t , < t < - , 
2 4 ' 



we find 



cos 3 1 



e(v) = (3m cost , e(v) — ve'iv) = (3m- 

cos 2t 



while from e'(v) = h(r) we get 



cos 2t 



msint 

Upon inserting these intermediate results in our principal envelope formula (2), we arrive at 
the "Martin-Roy-based" Coulomb lower bound 



E > min 

0<t<7r/4 



„ cos 3 1 T / cos 2t 
(3m ^ + V 



cos 2t 



msint 



5 Concave transform of harmonic-oscillator potential 



For the class of potentials V(r) which are concave transforms V(r) = g(h(r)), g" < 0, of the 
harmonic-oscillator potential h(r) = r 2 envelope upper bounds on the eigenvalues E n i of the 
spinless-Salpeter Hamiltonian (1) may be found for all energy levels (cf. Eq. (18) of Ref. 



E n£ < min 

r>0 



/ P 2 

(3\lm 2 + — + V(r) 



P = P n£ (2) = 2n + £-- 



(17) 



6 Concave transform of linear potential 

For the class of potentials V(r) which are concave transforms V(r) = g(h(r)), g" < 0, of the 
linear potential h(r) = r the "harmonic-oscillator based" envelope upper bounds (17) on the 
eigenvalues E n £ of the spinless-Salpeter Hamiltonian (1) may be (significantly) improved by 
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considering the composition of the envelope approximations for the kinetic-energy operator 
and the interaction-potential operator along the following lines. In the present investigation 
we focus our attention to attractive potentials V(r), that is, to all potentials which exhibit a 
monotone increase with increasing radial coordinate r: g' > 0. Every function g(r) that is a 
monotone increasing [g'(r) > 0] and concave [g"(r) < 0] function of r is a concave function 
f(r 2 ) of r 2 : g(r) = f(r 2 ) with f"(r 2 ) < 0. Consequently, all potentials V(r) studied in this 
section also belong to the class of potentials considered in Sec. 5. 

The square-root operator k(p 2 ) = \Jrri 1 + p 2 of the relativistic kinetic energy, regarded 
as a function of p 2 , is monotone increasing (i.e., k'i^p 2 ) > 0) and concave (i.e., k"(p 2 ) < 0). 
Let us assume that the interaction potential V(r) exhibits a similar behaviour, that is, that 
the operator V(r) = g(r) is monotone increasing, i.e., g'(r) > 0, and concave, i.e., g"{r) < 0. 
These operators are bounded from above by their respective upper-tangent approximations: 

HP 2 ) <ai(g) + 6i(?)p 2 , 



with 



g(r) < a 2 (t) + b 2 (t)r , 
bi(q) = k'(q) , a^q) = k(q) - qk'(q) , 

b 2 (t)=g'(t), a 2 (t)=g(t)-tg'(t) , (18) 

where p 2 = q denotes the point of contact between the kinetic-energy operator k(p 2 ) and its 
upper-tangent approximation a\(q) + bi(q)p 2 , and r — t labels the point of contact between 
the interaction-potential operator g(r) and its upper-tangent approximation a 2 (t) + b 2 (t)r. 
Consequently, our semirelativistic Hamiltonian H in Eq. (1) satisfies the operator inequality 

H = Pk{p 2 ) + V(r) < H(q, t) = (3 [ ai (q) + h^p 2 } + a 2 (t) + b 2 (t)r . 

The eigenvalues E n £(q, t) of the upper "tangential Hamiltonian" H(q, t) are given in terms of 
the numbers P n t{l) listed in Table 1 by the exact formula 

1/3 



E nt (q,t) = 3 



(3h{q)bl{t)P 2 t {l) 



+ p ai (q) + a 2 (t) . 



(19) 



In order to find the best such upper bound, we must minimize with respect to both q and t. 
The critical equations immediately yield the necessary conditions 

i 9'(i) q" 2 P 2 e (l) 

2(3k'(q) P ne (l) P ' 1 ) 

From the second of the above equalities we deduce that, for given quantum numbers n and £, 
the minimizing values of the parameters q and t are constrained to the curve qi 2 = P 2 ^!) in 
the (q, t) parameter plane. On this curve, and with the coefficients (18) of our upper-tangent 
approximations and the necessary (critical) conditions (20), the minimum of the expression 
on the right-hand side of the upper energy bound (19) may be cast into a new, simpler form: 




By remembering the above definitions of k(p ) and g(r), this expression eventually becomes 

<min .}\j m *+-— + Y(r) ■ P = P nl {l) . (21) 



E„ 



These upper bounds on the entire discrete spectrum {E n i\ are valid for any central potential 
V(r) which is a monotone increasing and concave function of the linear potential h(r) = r. 
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7 Variational upper bounds 

The standard tool for the derivation of rigorous upper bounds on the eigenvalues E of some 
self-adjoint operator H is the Rayleigh-Ritz variational technique [10, 11], which is based on 
the minimum-maximum principle: If the (discrete) eigenvalues E k , k — 0, 1, 2, . . . , of H are 
ordered according to E < E\ < E 2 < • • • , then the first d of them are bounded from above 
by the d eigenvalues Ek, k — 0, 1, . . . , d— 1, ordered according to Eq < E\ < ■ ■ ■ < Ea-i, of 
that operator H which is obtained by the restriction of if to a (/-dimensional subspace Dd of 
the domain of H, that is, E k < E k for all k — 0, 1, . . . , d — 1. Hence, for the accuracy of the 
obtained results, an appropriate, "reasonable" definition of this trial subspace Dd is crucial. 
If this d-dimensional subspace D d is spanned by a chosen set of d linearly independent basis 
vectors \ipk), k — 0,1, ... , d—1, the set of eigenvalues E may immediately be determined by 
diagonalizing the dxd matrix ({ipi\H\ipj)), i, j = 0, 1, . . . , d—1, that is, as the d roots of the 
characteristic equation 

det((^|#|^> - E = , i,j — 0,1, ... ,d — 1 . 

The relativistic virial theorem derived in Ref. [18] — for a very comprehensive review, see 
Ref. [19] — allows to define a precise quantitative measure [8, 9] for the quality of the results 
found within the framework of variational techniques: according to the analysis presented in 
Refs. [8, 9], for some generic Hamiltonian operator H consisting of a momentum-dependent 
kinetic-energy term T(p) and a coordinate-dependent interaction-potential term V(x), i.e., 

H = T(p) + V(x) , 

the accuracy of trial states \ip) which approximate the exact bound state under study can be 
quantitatively estimated by the deviation from zero of the quantity (called v in Refs. [8, 9]) 

(y|p-j§r(p)|y) t _ 
Q ~ <y,|x.JM/(xM { ~ V) - 

The main advantage of this measure for the accuracy of approximate eigenstates \tp) is that 
it does not require any information on the solutions of the investigated eigenvalue problem 
other than the one provided by the variational approximation technique itself. 

In order to get a first idea of the position of the lowest energy level E , a one-dimensional 
trial space is sufficient; a basis for the latter might be any of the two-parameter (normalized) 
trial functions 



0(x) = 



no?l n ( a 
■ exp — — r 



Here the scale parameter a > is varied. The parameter n > allows for later optimization. 
For the sake of definiteness, let us consider an interaction potential V(r) which is the sum of 
(attractive) pure power-law terms a(q) sgn(g)r 9 with q ^ and a logarithmic term a(0) lnr: 

V(r) = a(q) sgn(g)r 9 + a(0) lnr . 

We assume, of course, that the potential coefficients a(q) > are not all zero. With the help 
of Jensen's inequality [20], and by a tricky redefinition of the scale variable a, the variational 
upper bound on the ground-state eigenvalue Eq of the Hamiltonian (1), resulting from this 
choice of 0(x) as basis vector of our one- dimensional trial space, may be cast into the form 



En < min 

— r>0 



pJm 2 + — + a (<?) sgn(q)(V(n, q)r) q + a(0) \n(V(n, 0)r) 



T ' 2 9^0 



(22) 
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with the abbreviations 



p <^> - 2 Kin IwJ ' q * ■ 



where ^ denotes the digamma function 



1 dr(z) 
rr*) d^ 



8 The "Laguerre" trial space 

Unless the eigenstate corresponding to some eigenvalue is already an element of the trial 
space Dd-, any upper bounds on the eigenvalues of an operator bounded from below may be 
improved by enlarging Dd to higher d, or by spanning Dd by more sophisticated basis states. 

For the class of spherically symmetric potentials V(r) studied here, a very popular choice 
for the basis states which span the rf-dimensional trial space Dd of the variational technique 
are "Laguerre" trial states, given in their configuration-space representation by [21, 22, 6, 7] 



^ Mm(x) = \ T^T+TpTk'+l/ ^ 1 ex P(-^)^f +2p) (2^)^ m (fi x ) • (23) 

Here L ( ^\x) are the generalized Laguerre polynomials (for the parameter 7) [15] , defined by 
the power series 

4*(.) = D-i)'(*t?)£ 

and orthonormalized, with the weight function x 7 exp(— x), according to 


0^ m (f2) are the spherical harmonics for angular momentum £ and projection to, depending 
on the solid angle Q and orthonormalized according to 

The trial functions (23) involve two variational parameters, p (with the dimension of mass) 
and p (dimensionless), which, by the requirement of normalizability of these functions, are 
subject to the constraints p > and 2p > —1. 

One of the main advantages of the choice (23) for the variational trial states is the easy 
availability of an analytic expression for their momentum-space representation, obtained by 
Fourier transformation of ^^(x): 



^fc,<?m(p) 



(2p)^P+W. (-i) e \p\ e £ (-1)* ( k + 21 + 2p 



T(2£ + p + t + 2)(2pY ( 2t + p + t + 2 p + t 3 p 2 



p 2 + At 2)(2^+t+2)/2 ^ 2 ' 2 ' 2' P 2 +A* 2 

X ^m(fip) , 



12 



with the hypergeometric series F, defined by 

T(w) ™T(u + n)T(v + n)z n 
F(u,v:w;z) = — — — . > — : . 

The choice of the Laguerre states as the basis of the trial space D d allows to calculate the 
matrix elements of (any linear combination of) pure power-law potentials V(r) = a(q)r q and 
under certain circumstances (for instance, for \x = m for either radial or orbital excitations) 
the matrix elements of the kinetic term on purely analytical grounds, and to arrive therefore 
at analytical expressions for the dxd Hamiltonian matrix ({ipi\H\ipj)),i,j = 0, 1, . . . , d — 1. 
The explicit algebraic expressions of these matrix elements may be found in Refs. [22, 6, 7]. 
Up to and including a trial-space dimension d — 4, the Hamiltonian matrix ((ip^H^j)) may 
be diagonalized algebraically. For the funnel potential (14), the case d — 1 yields the bound 

'64/3 
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9 The "local-energy" theorem 

Some information on both upper and lower bounds on the isolated eigenvalues E n £ of some 
Hamiltonian H may be gained with the help of the so-called "local-energy" theorem [23, 11]. 
Unfortunately, the proof of this fundamental criterion makes use of the "nodal theorem" for 
the eigenstates of H. As a consequence of this, for Hamiltonians H defined in more than one 
dimension the local-energy theorem can be applied only to the ground state of the spectrum. 
It has been applied to the (three-dimensional) "relativistic Coulomb problem" in Ref . [24] . 

Upper bounds on the energy levels of Hamiltonians bounded from below can be obtained 
with considerably more efficiency by the variational technique described in Sec. 7. Therefore 
we employ here the local-energy theorem for the derivation of lower bounds on the spectrum 
of H. In order to formulate the local-energy theorem in momentum-space representation, we 
introduce the "local energy" 



, d 3 gU(p-q)0(q) 

S(p) = pJm 2 + p 2 + 

Here, V(p) is the Fourier transform of the interaction potential V(x) under consideration: 

^ (p) = jhf I d * x • x ) v w • 

</>(p) denotes a suitably chosen, positive trial function: 0(p) > 0. The "lower-bound part" of 
the local-energy theorem then states that the lowest-lying eigenvalue E of the Hamiltonian 
H is bounded from below by the infimum of £(p) with respect to the momentum variable p: 

£ >inf£(p). (24) 

The application of the local-energy theorem to confining potentials demands special care 
for the following reason: the Fourier transform of any confining potential is not well-defined; 
the potential has to be regarded as a distribution and an appropriate regularization must be 
applied. For instance, the linear potential Vh(r) = Xr, where A is a coupling parameter with 
mass dimension 2, may be regularized by the exponential exp(— er), e > 0, or by writing it 
as the derivative of this exponential or as the second derivative of a Yukawa-type function: 

d 3^ exp ( — <st) 

Vh(r) = Xr = Alimrexp(— er) = — Alim — exp(— er) = Alim 



eio " " .0 de £|0 i); 



2 
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With these regularizations, the Fourier transform Vl(p) of the linear potential Vl(t-) reads 

r,.!p) r,.ip) 4ii"' I! " _ ?y: ~ 



7T 2 elO (p 2 + e 2 ) 3 

Note that (as a consistency check) the integral of the Fourier transform VL(p) has to vanish: 
J d 3 p Vl(p) = 0. This holds for every power-law potential V(r) = a n r n with exponent n > 0. 
The singularity structure of VL(p) becomes manifest when writing the derivative in Vi,(r) as 
a difference quotient, 

Vj(r)=Ar = -Alim eXp( - er) - 1 , 

ej.0 e 

which yields 

"1 



^L(p) = VL(|p|) = Alim 

ej.0 



;^ (3) (P)- 



e 7T 2 (p 2 +£ 2 ) 2 

Thus the Fourier transform VL(p) of the linear potential consists of a (negative) regular part 
and a singular part, which is a distribution localized at the origin p = 0. However, a detailed 
inspection shows that for the case of the funnel potential the momentum-space ground-state 
eigenfunction is indeed positive. This fact justifies our use of a positive trial function 0(p). 
The simplest normalizable trial function which comes to one's mind is the exponential: 

0( r ) = ^- ex P(~^ r ) • 

This choice for the ground-state trial function <fi(r) corresponds to the case k = £ = m = 
and p = 1 of the Laguerre basis functions ipk^mi^) introduced in Sec. 8. It is straightforward 
to calculate the Fourier transform of 0(r), or to extract it from the general expression given 
explicitly in Sec. 8: 

,/ \ ,„ ,v . ( n |p|\ 2v/27? 

P = P = | [7 2 , — 2T sm 2 arctan — = — — — — . 

7r|p|(p 2 + /i 2 ) V /I J 7r(p 2 + /i 2 ) 2 

In order to illustrate the application of the local-energy theorem, let us consider the example 
of a potential V(r) which is a sum of pure power-law terms: 

V(r) = J2 a (<l) rq ■ 

i 

For this power-law potential and the ground-state trial function 0(r), it is straightforward to 
obtain the corresponding local energy £(p) by Fourier transformation of the product r 9 0(r): 



[q + 2) arctan 



IPI 



The variational parameter \i entering into our ground-state trial function 0(r) either may be 
given a chosen fixed value, or it may be adjusted in order to maximize the lower bound (24). 



10 Sums of distinct potential terms 

For the class of potentials V(r) which consist of two or more distinct terms V^(r) ^ V^\r) 
for % ^ j, where, for every single component problem defined by the "one-term" Hamiltonian 
Hi = P^ni 2 + p 2 +V^(r), information about the bottom of its spectrum is available, these 
individual pieces of information may be combined to form a lower bound on the spectrum of 
the Hamiltonian (1) with such a sum potential V(r). The precise requirements for this "sum 
approximation" to work are the following: 
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1. The interaction potential V(r) entering in the semirelativistic Hamiltonian of Eq. (1) 
is a sum of (more than one) different components V^\r) = CihS l \r): 

i 

2. Every component problem Hi = (5\/m 2 + p 2 + q/iW (r) supports, for sufficiently large 
values of its coupling parameter q, a discrete eigenvalue at the bottom of its spectrum. 

The sum approximation has already been studied for the case of Schrodinger Hamiltonians 
(involving nonrelativistic kinetic energies) in Refs. [1, 25]. The analysis of the corresponding 
semirelativistic problem has been presented in full generality in Ref. [26]. Therefore, here we 
summarize only the result for the case in which the potential V(r) is a sum of attractive pure 
power-law terms: 

V(r) = 5>(g)sgn(g)r«, (25) 

where the coupling parameters a(q) are non-negative and not all zero. For this potential, a 
lower bound to the bottom of the spectrum of the semirelativistic Hamiltonian (1), which is, 
by assumption, the lowest eigenvalue E of H, is given by the expression 



En > min 

~~ r>0 




(26) 



+ XX?) sgn(g) (P(g)r) 

if, for the particular potential V(r) under consideration, there exist (or there can be found) 
P numbers P(q) such that, whenever there is only one term present in the potential (25), the 
right-hand side of the inequality (26) yields either the exact ground-state energy eigenvalue 
of the corresponding "single-term" problem or a lower bound to it. 

Consequently, when applying our general sum lower bound (26) to a given sum potential 
V(r), the main task is to derive appropriate P numbers for all components in the sum (25). 
Let us illustrate this by considering the two examples relevant for the funnel potential (14). 

• For a Coulomb component potential h^(r) = — 1/r, corresponding to the case q — — 1 
in Eq. (25), we may take advantage from the explicit results presented in Sec. 4. By the 
change of variables r — > Pr, the Coulomb lower bound (12) may be cast into the form 



En > min 

~~ r>0 




, v < (3v P . 



V{Pr) 

Specifying this general result to the Coulomb potential V(r) = —v/r yields the bound 
Eq > min 



r>0 



v < (3v P 



This bound is, of course, nothing else but, for P = 2/ir, the Herbst lower bound (9) or, 
for arbitrary P, the Herbst-like lower bound (11) to the relativistic Coulomb problem. 
Consequently, for the Coulomb component term we may identify the P number P(— 1) 
required for the lower bound (26) with the lower-bound parameter P defined in Sec. 4: 
P(-l)=P 

For a linear component potential h®(r) = r, that is, for the case q = 1 in Eq. (25), (a 
possible value of) the P number P(l) required for the lower bound (26) turns out to be 
related to the ground-state eigenvalue e (v) of the Hamiltonian operator ^/p 2 + vr by 

e (v) = 2jvP(l) ; 

a numerical determination of this lowest energy eigenvalue en{v) yields P(l) = 1.2457. 
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11 Comparison: The case of the funnel potential 

In order to demonstrate the power of the different energy bounds derived so far, let us apply 
the above results to the very illustrative example of the Coulomb-plus-linear potential (14). 
By factorizing off an overall coupling constant < v < 1, we write this potential in the form 

V(r) = v(-^ + br^j . (27) 

For the coupling parameters entering in this funnel potential V(r), different notations have 
been used in Eqs. (14), (25), and (27). Evidently, the three sets of potential coefficients have 
to be identified according to a(— 1) = C\ = av > and a(l) = c 2 = bv > 0. 

Fig. 1 compares the various energy bounds derived or investigated in Sees. 4 through 10 
for the ground-state energy eigenvalue Eq of the Hamiltonian (1) with the funnel potential 
(27) as a function of the overall coupling parameter v. These energy bounds necessarily form 
a bunch of concave curves E(v), all starting at the (free-energy) value Eq(v = 0) = f3m = 1. 

2.1 
2 
1.9 
1.8 
1.7 
1.6 

E(v) 

1.5 
1.4 
1.3 
1.2 
1.1 
1 

0.2 0.4 0.6 0.8 1 

v 

Figure 1: Comparison of bounds to the ground-state [(n, £) = (1, 0)] energy eigenvalue of the 
semirelativistic Hamiltonian H = j3y/m 2 + p 2 + V(r) with a Coulomb-plus-linear potential 
V(r) = v(—a/r + br), for a = 0.2, b = 0.5, m — f3 — 1. The bounds are plotted as functions 
E{y) of the overall coupling parameter v in V(r). The energy bounds compared here include 
(from top to bottom) the harmonic-oscillator upper bound (H) of Eq. (17), the linear upper 
bound (L) of Eq. (21), the "one-dimensional" variational upper bound (V) of Eq. (22) with 
n = 1.74, the variational upper bound (E) based on a 25-dimensional Laguerre trial space, 
which should come rather close to the exact eigenvalue, the sum lower bound (S) of Eq. (26) 
with P numbers P(— 1) = 0.72811 and .P(l) = 1.2457, the lower bound (T) provided by the 
local-energy theorem (24), and two Coulomb lower bounds (C) according to Eq. (12), where 
with respect to the parameter P the lower curve corresponds to the fixed value P = 0.72811 
whereas the upper curve corresponds to the optimized value of P calculated from Eq. (16). 
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The Coulomb lower bound (12) may be applied to the funnel potential either with a fixed 
value of the lower-bound parameter P or with the optimized P(v) computed from Eq. (16). 
For "fixed-P" bounds, the maximum value of P allowed by the Coulomb coupling constant 
constraint discussed in Sec. 4 is found, for given values of the kinetic-term parameters ft, m 
and potential coefficients Ci, c 2 , as solution of that relation which is obtained when equating 
the left-hand and right-hand sides of the inequality (15) with the critical coupling vp given 
by the definition (13), that is, by solving, for the maximum values of c\ and c 2 occuring here, 
the constraint (16): 

ci + j^—^ppVT^. 

1 — F z m z 

For a single, unit-mass particle, i.e., (3 — 1 and m — 1, and the values of the funnel-potential 
coupling parameters a and b used in Fig. 1, a — 0.2 and b = 0.5, this yields P = 0.728112. 

We try to approach the state we are interested in as closely as possible by the application 
of the Rayleigh-Ritz variational technique briefly sketched in Sec. 7, with the trial space Dd 
spanned by the Laguerre basis defined in Sec. 8. By fixing the variational parameter p to the 
value p = m, we are able to obtain the matrix elements of the Hamiltonian (1) analytically, 
avoiding thereby the necessity to use a numerical integration procedure for the evaluation of 
these matrix elements; for definiteness, we fix the variational parameter p to the value p — 1. 

The accuracy of such an approximation to an eigenstate may be estimated with the help 
of the "virial-theorem-inspired" quality measure Q discussed in Sec. 7. Fig. 2 shows, for the 
trial-space dimension d = 25 chosen for the variational bound (E) in Fig. 1, the dependence 
of the accuracy Q on the overall coupling v. We obtain Q — \Q\ < 1.1 x 10~ 4 over the whole 
parameter range < v < 1. Such high accuracy should be sufficient for the present purpose. 
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Figure 2: Measure Q of the accuracy of the (variational) approximation of the ground-state 
eigenfunction of the semirelativistic Hamiltonian of Fig. 1 [that is, H = f3\/m 2 + p 2 + V(r), 
where V(r) is the Coulomb-plus-linear potential V(r) = v(—a/r + br) with a = 0.2, b = 0.5, 
m = (3 = 1] by a superposition of 25 Laguerre basis functions (23) with p — m and p — 1. 
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For the sake of completeness, we plot, in Fig. 3, the behaviour of the ground state of the 
Hamiltonian (1) with the Coulomb-plus-linear potential (27), for the maximum value v — 1 
of the overall coupling parameter v, determined by the variational technique of Sec. 7, with 
the Laguerre basis of Sec. 8, in both configuration and momentum space. As expected, the 
ground-state momentum-space eigenfunction possesses no node, which is in agreement with 
our assumption of a positive momentum-space trial function 0(p) (which is also depicted, as 
the case d — 1, in Fig. 3) in the derivation of the local-energy lower energy bound in Sec. 9. 




Psi 




Figure 3: Variational approximation of the lowest-lying eigenstate in configuration space (a) 
and momentum space (b) of the semirelativistic Hamiltonian (1), H = fi^/m 1 + p 2 + V(r), 
with a Coulomb-plus-linear potential V(r) = —a/r+br, for a = 0.2, b = 0.5 and m — (3 — 1, 
by a superposition of Laguerre basis functions (23) with the variational parameters fixed to 
/i = m and p = 1, for the trial-space dimensions d — 25 (full lines) and d — 1 (dashed lines). 
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12 Summary and conclusions 



This review of the discrete spectra of semirelativistic Hamiltonians H of the form (1) should 
be considered as a compilation of rigorous (semi-) analytical results for both upper and lower 
bounds on the discrete eigenvalues E n g of H, presented in a way suitable for immediate use: 
the algebraic structure and the convexity properties with respect to suitable basis potentials 
h(r) of the interaction potential V(r) one wishes to use in one's semirelativistic Hamiltonian 
(1) unambiguously determine which of all the energy bounds discussed here actually apply. 
With the explicit analytical expressions for these bounds on the discrete eigenvalues at one's 
disposal, one is left with the (numerical) optimization with respect to a single real variable. 
Clearly, within the respective limitations, the Rayleigh-Ritz variational technique of Sees. 7 
and 8 as well as the local-energy theorem of Sec. 9 apply to arbitrary reasonable potentials. 

For the particular example considered in Sec. 11, viz., the Coulomb-plus-linear potential 
(27), the sum lower bound turns out to be somewhat better than the lower bound provided 
by the local-energy theorem. Clearly, the sum lower bound cannot always be better than the 
local-energy bound because one could, in principle, choose in the local-energy theorem (24) 
the exact eigenfunction for the trial function <fi and thus recover the exact energy eigenvalue. 

In conclusion, let us stress that all but two of the various energy bounds we have studied 
in this work may eventually be expressed by essentially the same general semi-classical form. 
We can view this in the following way: Since the square of the momentum p 2 scales like 1/r 2 , 
we expect that 




where P is a suitable coefficient; thus the balance between the kinetic and potential energies 
required by the minimum-maximum principle might lead to a search for the minimum of an 
energy expression of the form 



What makes this heuristic argument interesting is the fact that definite values of P can be 
prescribed so that this minimum provides a priori a bound on the (unknown) exact energy. 
In the case where V(r) is a sum of terms, a trivial change of variables is used to move the P 
parameters inside the potential V(r), and then optimal distinct parameters can be used for 
each potential term. 

The lower bound by the sum approximation incorporates our best known lower bounds 
for the sub-problems containing only a single potential term; the final result is the optimized 
mixture of these contributions. The lower bound by the local-energy theorem (24) could, in 
principle, achieve arbitrarily high accuracy. The difficulty here, as with variational methods, 
is to find a suitable trial function. To the extent that we have searched we did not manage to 
find a better lower bound with the help of the local-energy theorem (24) than that provided 
immediately by the sum approximation (26). 

Our most accurate results overall were upper bounds obtained by using a 25-dimensional 
trial space spanned by Laguerre functions. We estimate the error for these results to be less 
than 0.01% over the total range of potential parameters studied; such accurate estimates are 
invaluable for estimating the effectiveness of the large variety of alternative simpler energy 
bounds which we have reviewed in this article. 
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